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I. INTRODUCTION 



Marcin Kazmierczak* 

Institute of Theoretical Physics, University of Warsaw ul. Hoza 69, 00-681 Warszawa, Poland 

In this article the methods of canonical analysis and quantization that were reviewed in the first 
part of the series are applied to the case of the Dirac field in the presence of electromagnetic inter- 
action. It is shown that the quantization of electrodynamics, which begins with a given Lagrangian 
and ends by perturbative calculation of scattering probability amplitudes, can be performed in the 
way that does not employ Poincare symmetry of space-time at any stage. Also, the causal structure 
is not needed. 
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' In this second part of the sequel, I shall apply the canonical formalism reviewed in the first part to the case of 
electrodynamics with spinorial matter. I shall begin with the discussion of constraints and gauge transformations in 
Section II. The issue of integrating infinitesimal gauge transformations to the finite form will be addressed and the 
^ c"| total and extended Hamiltonian formalisms confronted. 

O |. In Section III, the issue of consistent imposition of gauge conditions will be discussed on the simplest example of 
Coulomb gauge. The generalized Dirac bracket (GDB) will be constructed for this case. Then equations of motion 
will be discussed. Since they will appear to be by far too complicated to be exactly solved, the transition to the 
interaction picture will be necessary. When this is done, the equations simplify immensely and straightforward Fock 
quantization can be applied to the interaction picture fields. Knowledge of this interaction picture representation will 
appear to be sufficient for the perturbative description of scattering processes. 

In Section IV, the two examples of Compton scattering and electron, positron — > muon, anti-muon scattering will 
. be discussed. The elements of the S matrix in the lowest non-trivial order (i.e. the second order in fine structure 
constant) will be explicitly computed. No Feynman rules will be postulated, although the relation of the calculations 
\Q . to the standard quantization based on Feynman diagrams will be explained. 

o\ 

(N 

II. THE STRUCTURE OF CONSTRAINTS OF ELECTRODYNAMICS 



The theory of electromagnetic interaction is interpreted as a gauge theory of a U{1) group. The Lagrangian for 
fermions in the presence of this interaction is obtained through the minimal coupling procedure. Explicitly, 

C = I ^Vai, - P^ 7 » - mW> - -F ab F ab = C D + C EM - eAjp^, 

><: \ _ _ 4 (h i) 

£d = - (ii>l a d a i> - d a *PY l iP) , C EM = -^F ab F ab , F ab = d a A b ~ d b A a , , V a i\> := (d a + ieA a )^, 

where e is the electric charge and m represents the particle's mass. In order for the covariant derivative T> a ip to 
transform as T> a 4> —> T>' a i\)' = e~ leX T> a 4> under the U(l) gauge transformation ip tjj' = e~' leX ijj, the U(l) connection 
one-form (electromagnetic four-potential) needs to transform as A a — > A' a = A a + d a X. Clearly, the Lagrangian is 
then invariant under these transformations. The fields of the theory are tp, tp and A a . The equations for conjugated 
momenta are 

<>C %°, W =-^=- = -U°^ «° = lJ^ = -F 0a , (H.2) 



d(doi>) 2" ' ' d{d Q i>) 2 ' r ' d(d A a ) 

from which the primary constraints follow 



Xl = 7T - 2^7°, X2 = TT + ^1°^P, ll=TT°- (IL3) 



'Electronic address: marcin.kazmierczakafuw.edu.pl 



2 



There is a possible source of confusion in the notation, since the latter 7 is used to denote the Dirac matrices, as well as 
first class constraints (I will argue that "fi is first class later on) . The confusion can be avoided if one remembers that 
the Dirac matrices are denoted by 7 with superscripts and the constraints are labeled by subscripts. The canonical 
Hamiltonian is given by 

H = j Hd 3 x, 

H = nip + 1jk + tt% - C (II.4) 
= ~4j (-ij j dj + m) ip + ^ttV + ^FijFij + eAjip^ip + A (etp^ip - djir 3 ) + d (a -k 3 + ^py>ipj 



and the total Hamiltonian is 

H T = H + 



J ( U (f) 7l (f)+ X i(x) M 1 (x)+u 2 (f) X2 (x)),d 3 a ; (II.5) 



where u's are arbitrary multipliers (note that u 1 and u 2 are matrices). The last component of H in (II. 4) is a three- 
divergence that results in a surface term in H . Such terms do not contribute to the functional derivatives of H with 
respect to fields, since it is common to consider variations that vanish on the boundary of the integration region when 
defining variational derivatives. If the fields themselves (not only their variations) vanish at spatial infinity, which is 
usually assumed in the case of flat space-times, then the surface terms simply vanish thus not giving any contribution 
to the total energy. 
(III.6) 

The GPB is given by 

r ( SF SG SF SG SF SG SG SF 5G 5F SF SG \ 3 

[ ' lGP ~J \SA a (x) Sir a (x) ~ 8^{x) SA a (x) + S^(x) Sir(x) Stt{x) 8^{x) T 5iP{x) dn(x) _ Stt(x) S$(x) J *' 

(11.6) 

where the upper sign applies whenever at least one of the variables F, G is even and the lower one corresponds to F 
and G odd (see formula (III. 6) of [14] and the discussion therein). The consistency conditions for the time evolution 
of constraints are 

[71, H T ] G p = ~y^t = ~yt~ = djir 3 - eipj°ip w 0, 
SA dA 

[Xi,H t ]gp = - ^^pV = -iVpPl - m^-eAo^ - iu 2 f « 0, (II.7) 

[X2, H t ]gp = ~^=- + ^7°^" = H 3 Vjip -imp- eA "f°ip + i^u 1 w 0. 

Olp 2 07T 

The first equation does not depend on u's and hence gives rise to a secondary constraint 

72 := djir 3 — eip~/°ip. (II.8) 
The second and the third equation of (II. 7) yield merely the restrictions on u's of the form 

U 1 = (-7°7 J 'X> J - imj ) ip - ieA ip, U 2 =U T = -Vftp^-y + im^~f° + ieAtfp. (II.9) 
The bracket of 72 with Ht is 

[72, H t ]gp = -e [dj (^7» + V^V + AV] ■ (H.10) 

This is yet another restriction on u's, but it is easy to verify that it is satisfied automatically if u 1 — U 1 and u 2 = U 2 
(see (II. 9)). Hence, neither farther restrictions on u's nor additional constraints are produced. The consistency 
algorithm is accomplished. 

Let us now investigate what class do the constraints belong to. Clearly 7 1 commutes with all the others and hence 
is first class. The brackets of the remaining constraints are nontrivial 

[12{x),xuS']gp = -e(ip(x)'y ) l 5(x-x'), [72(2), X21x>]gp = -e (j ip{x)) ^(x ~ x"), 
[Xiis, X21'S']gp = i^i n S(x-x / ). 
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It seems that there are three families of second class constraints, xi> X2, 7 and one family of first class ones, 71. 
However, as explained at the end of Sec. II. A of [14], the constraints are well separated if the number of first class 
constraints is possibly large, i.e. no first class constraints are hidden in the second class ones. Are these constraints 
well separated? To verify this, let us try to construct a first class constraint 7 from a linear combination of Xi's, X2' s 
and 7's 
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:= J ( X i(x)Ai(f) + A 2 (x)x2(x) + K (f)7 2 (f))d 3 x, (11.12) 



where A's are matrices and n a scalar function. In order for this constraint to commute with all the others, it is 
necessary the the following relations between Ai, A2 and n hold 

i/>7°Ai « — A27 < V, A2 ~ —ieip, Ai s=s ienip. (11.13) 

They are sufficient for 7 to be first class. Hence, a first class constraint 

7 = y k(x) (72 (x) + iexi(x)ip(x) - ieip(x)x2(x)) d 3 x (II. 14) 

can be constructed. Since n remained completely arbitrary, the family of first class constraints 

72O?) := 72O?) + iexi(x)i)(x) - ieip(x)x2{x) (11.15) 

was constructed. Since the system of constraints xi> X2, 71, 72 is equivalent to xi, X2, 71, 72, we are allowed to use 
the latter. To summarize, the constraints of the theory are 

7i=7T°, 72 = d-jn 3 + ie(mp - ipT), Xi = ^ ~ , X2 = 7f+^7°V', (11.16) 

where 7's are first class and x's are second class. These constraints are well separated, since it is impossible to 
construct a first class constraint from linear combination of Xi's and X2' s - Using (II. 9), the first class Hamiltonian 
H 1 can be calculated 



H' = H + J (xiU 1 + U 2 X2 ) d 3 x 



ip (-i-fWj + m) tp + ^TT j n J + ^FijFij - A j 2 + Xi {-"p^V^ - im-f°ip) + (-Vjip^j + irmp'y ) X2 
The total and extended Hamiltonians are then 



d 3 x. 
(11.17) 



H T = H' + J w (x) 7l (f) ( i 3 x, H E =H' + J (u°(f)7 1 (x)+w;(x)7 2 (f))d 3 x, (11.18) 
where u° and w are arbitrary functions. 
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A. Equations of motion 

The equation of motion of any dynamical variable F generated by He is F = [F,He]gp- For basic canonical 
variables this gives 



Aq = Uq, 


(//.19a) 


Aj = tt 3 +dj(A -w), 


(77.196) 


7T° = djW-i + ie(Trip — ipw) = 72, 


(//.19c) 


* 7 n 1 " / 7 / 1 ~ T 7 — \ 

n J = Oi-tij + ie(TTj Yyj + yjj 7 j tt), 


(11 .19a) 


t/> = (— 7°7 J 'fj — 17717°) "0 — *e(Ao — w)^; 


(77.19e) 


t/> = —Vjip-fjj + 77771/17° + 7e(^4 — w)ip = i>, 


(77.19/) 


7T = —Dj'K^j ^ + imirj + 7e(^4o — w)n, 


(77.195) 


TT = (— ^ ^ 3 T>j — 77777°) 7f — * e (^0 — = 71"; 


(77.19/i) 


7T°=0, 


(77.197) 


djir 3 = ie(ipn — Trip), 


(7/.19j) 


n = ^7°, 


(//.19fc) 


i , 
t = -g 7 ^' 


(7/.19Z) 



(11.19) 



Vjtp := (9j + ieAj)ip, Vjtt :— (dj + ieAj)-K, Vjip := Vjip, Ujir :=Vjir. 



The first four equations a,b,c,d,e, f,g,h follow from evaluation of GP of the variables with He and the equations 
k, I are the constraints. A high degree of redundancy can be seen in this system. The equations i, k, I and 6 can be 
used to eliminate 7r°, tt, tt and tt j from the remaining ones. Then (II. 19a) just states that Ao is an arbitrary function, 
(II. 19c) becomes equivalent to (II.19j), (II.19d) gives 

8 a F aj = e^ip + djw, (11.20) 

(II.19e) gives 

(ij a V a -m)ip= -ewj°iP, (11.21) 

(II. 19f ) follows from (II.19e) by Dirac conjugation, (II.19g) and (II.19h) reduce to (II.19f) and (II.19e) on the constraint 
surface. Finally, after elimination of tt, tt and tt- 7 , (II.19j) gives 

d a F a0 = e^yV - didiW. (11.22) 

In the case when w — 0, (11.20) and (11.22) are Maxwell equations and (11.22) reduce to the £/(l)-co variant Dirac 
equation of spinor electrodynamics. In the w ^ case the equations, when written in terms of A a , do not assume the 
usual form of Maxwell equations. This example reflects a general rule that it is the dynamics generated by Ht, and 
not He, which is equivalent to the Euler-Lagrange equations of the initial Lagrangian (see the discussion below the 
formula (11.20) of [14]) . However, if A is transformed into A := A — w then, in terms of the quantities F a b and 
Vtp that contain A in place of A , the equations assume the standard form 

d a F ab = e^'ip, (i"f a V a - m) t/> = 0. (11.23) 

Whether these equations are equivalent to Maxwell equations or not depends on how the physically measurable electric 
and magnetic fields are defined in terms of A a ■ The correct definition in the extended formalism is 

E l ■= tt 1 = F M + d iW = F ia , B { ■= e ijk djA k , (11.24) 

where e ljk is antisymmetric symbol with e 123 = 1. Note that the position of spatial indexes is important, since in the 
metric convention (+, — , — , — ), which is used, shifting a spatial index leads to the change of sign. 
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B. Gauge transformations 

The time evolution of any dynamical variable F is given by 

F=[F,H e }gp = [F,H') G p + J ( U °(f)[ J F 1 ,7i(a ? )]GP + ^(2 f )[^72(2 ? )] G p)rf 3 a ; (11.25) 

and depends on arbitrary functions u° and w. The transformations corresponding to changes in these functions ought 
to be interpreted as gauge transformations, as explained in Sec. II. B of [14]. From (11.25) it is clear that the action of 
a general such transformation on F is 

SF = J (Su°(x)[F, 7l (x)} GP + 5w{x)[F, l2 {x)] GP ) d 3 x. (11.26) 

The transformations corresponding to the change in u° are said to be generated by the constraint 71, whereas those 
following from the change in w are generated by 72. Imagine that F(t) and F(t) correspond to the dynamics obtained 
for the choice of arbitrary functions as u°, w and u° = u° + Su°, w — w + Sw respectively. Assume that F(to) — F(to) 
at some instant to- The difference SF(t) = F(t) — F(t) corresponds to the unphysical gauge transformation. But how 
can this difference be calculated explicitly? If t = to + r and F is analytic, then 

5F (t +r) = J2SF^(t )-, (11.27) 

n=l 

where F^> denotes the n-th time derivative of F. Up to first order in r, this change can be easily evaluated by means 
of (11.26). But what if we wished to find a finite form of a gauge transformation? 

The equations of motion in the form (11.19) are very helpful. First note that instead of using (11.26) it is straight- 
forward to obtain the variations of first derivatives of basic canonical fields directly from (11.19). From now on, we 
will be interested in changes in the derivatives at t = to but the argument to will not be written explicitly. So for 

example SF := F(t ) — F(to) = for F being any canonical variable, but SF := F(to) — F(to) will not vanish in 
general. Specifically, from (11.19) it follows that 

SA = 5u°, SAj — —djSw, 5-k a — S-k^ =0, 5ip = ieSwip, Sip = —ieSwtp, 5-k = —ieSw-K, 5n — ie5vm. 

(11.28) 

To obtain variations of higher derivatives, differentiate each side of (11.19). Then the second time derivatives will 
be expressed by the fields and their first derivatives, whose variations are already known. Variation of the resulting 
system leads to 

5Ao = u°, 8A 3 = d 3 {5u° - Sw), , = -ie(Su° - Sw)tP - e 2 (8w) 2 i) + 2ieSwip, , 5tt> = (11.29) 

(variations of the derivatives of the remaining fields can be found by conjugations and application of the constraints). 
Note that Su° and Sw are not infinitesimal and therefore some care is required in calculating the variations. One 

should simply calculate the difference, e.g. Sij) = ip — "ip, so that the terms proportional to the higher powers of Su° 
and Sw are not omitted. 

The variations of higher derivatives can be found by continuation of this iteration procedure. The Tylor series that 
emerges can by shrunk to 

5A (t,x) = d t X(t,x), 5A 3 {t,x) = dj (X(t,x) + n{t,x)) , 5ip(t,x) = (^-MHt,s)+ K (t,s)) _ ^ 

f t ft' ft (11.30) 

\(t,x):= dt' dt"Su°(t",x), n(t,x):=- dt'Sw(t',x). 

J to J to J to 

The transformation of ir^ can be found from (II. 19b). It follows that 8^{t,x) = 0. Hence, it is ir^ = F-? + djw, 
and not F^ , which is gauge invariant and hence can be interpreted as a measurable physical electric field EK If 
this interpretation is assumed than the equations of motion that follow from the extended formalism, when expressed 
in terms of measurable quantities E, B, j a — eip^ip are precisely the Maxwell equations (note that from (11.30) it 
follows that j a is indeed gauge invariant). 
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III. COULOMB GAUGE 



Now that we have all the gauge freedom explicitly described, we should proceed to quantize the theory. Since the 
only second class constraints are xi and X2> the generalized Dirac bracket is given simply by the formula (III. 23) of 
[14], although the generalized Poisson bracket is now given by (II. 6). If GDB's are promoted to the commutators of 
the operators in the quantum theory, the second class constraints can be consistently interpreted as strong operator 
equalities. Were there no first class constraints, we could try to quantize the theory in much the same way as we 
did for the free Dirac field in [14]. But the first class constraints are the obstacle. They cannot be interpreted 
as strong operator equations, since their GDB's with other dynamical variables do not vanish in general. One 
method of implementation of these constraints in the quantum theory is to demand that the physical Hilbert space of 
physical quantum states is a subset of a larger kinematical Hilbert space defined by the condition j a \il>) = 0, where 
a enumerates all the first class constraints. This method was originally proposed by Dirac [2] and gave rise to the 
BRST quantization. However, I shall use a different method, which allows the computation of physically measurable 
quantities most quickly. This is the fixed gauge quantization. In this approach, one simply adds another constraint 
to the theory, called gauge condition, which results in all the constraints being second class at the end. A gauge 
condition is any relation between the q's and the p's of the form 

X (q,p) « (III.l) 

which is accessible, i.e. any point (q,p) can be transformed by gauge transformation to the one that satisfies (III.l). 
Also, we wish that the gauge condition, after being inserted to the consistency algorithm together with other con- 
straints, eliminate all the gauge freedom. Otherwise, we need to impose farther gauge conditions or use the Dirac 
quantization anyway. 

A shall consider the simplest condition, which is the Coulomb gauge 

X = dtA 1 = VA« 0. (III.2) 

This is clearily accessible by the transformation (11.30), since for any A one can construct A' = A — V/, / = A + k 
that satisfies V 'A — 0. One only needs to take 

f(t,x) = -±- / ^ ( *'f[W (III.3) 
J y ' ' A-k J \x - x'\ K ' 



To see this, use the identity 



-A^-i- = V* ( |t ) = 4n5 3 (x - y). (Ill A) 



\x-y\ \\x-y\ 



Let us now step back to the point of the analysis in which the primary constraints were just found, eq. (II. 3). We 
shall now add the gauge condition (III. 2) to the set of primary constraints (the reader is encouraged to find out that 
the results would not alter if x was interpreted as a secondary constraint). Hence, the primary constraints and the 
total Hamiltonian are now 

Xi = 7T- ~^7°> X2=tT+^7V, 7i= A X = djA\ 

Z . 1 (III.5) 

H T = H+ / (u (f) 7l (f) + X i(x)u\x) + u 2 (x) X 2{x) + u{x)x{x)) ,d 3 x, 

where H is still given by (II. 4). Now that the set of primary constraints and Ht have changed, it is necessary to 
rerun the consistency algorithm for the time evolution of constraints. As before, the bracket of ir° with Ht gives the 
constraint 72 and requiring the brackets of xi and X2 with Ht to vanish leads to the restrictions on u 1 and u 2 (II. 9). 
However, calculation of the bracket of x with Ht results in a new constraint 

SHt 

[X,H T ] GP = [djA^H^GP = -dj-^T = ~ d -j* J - AA o 0. (III.6) 

So, in the first run of imposing consistency conditions the two secondary constraints 72 and (f> were found and the 
restrictions (II. 9) on u 1 and u 2 . Continuing with the algorithm for 72 and 4> one gets 



[72, H t ]gp = ~e [dj (W» + ^"w 1 + « 2 7>] + djdju w 0, 
[<j>, H t ]gp = edj (^ 7 V) + djdj(u - u°). 



(in.7) 
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If (II. 9) is used (which is allowed), the first equation gives the restriction on u 

Au = 0. (III.8) 

The second equation than gives 

Ait = edj (ipy'ip) . (III.9) 

Hence, no farther constraints were produced and the consistency algorithm is now finished. Instead of 72 the constraint 
72 = djiri + ie(irip — tpw) can be used, as before. The commutation relations of the constraints are now the following: 
7T° commutes with everything but (f> 

[ir°(x),cf>(y)} G p = A6(x-y). (111.10) 

72 commutes with everything but \ 

h2(x), X (y)]GP = -AS(x-y). (III.ll) 

Xi and X2 commute with all the other constraints and the bracket of xi with X2 is the same as before. Finally, the 
bracket of x with <f) is 

[x(x-),(t>(y)} G P = -A5(x-y). (111.12) 

Hence, all the constraints are now second class, as desired. It is straightforward to check that these constraints do not 
hide any first class ones, i.e. it is impossible to construct a first class constraint from a linear combination of them. 



A. The generalized Dirac bracket in the Coulomb gauge 

Before calculating the GDB, it is useful to rearrange slightly and rename the constraints. In the following I shall 
use 

xi = K--^°, X2 = tf+27V, 

X3 = 72 =d 3 7T : > +ie(irip -ipT), x± = X = djA 3 , (111.13) 
X5 = 71=^°, xe = <f> + 72 = ie(Trip - ijm) - AA . 

Hence, instead of <f> I will now use x& — 4> + 72- This is because now the constraints can be grouped into pairs Xii Xi\ 
X3jX4; X57X6; such that the constraints in a given pair have vanishing brackets with those in the remaining pairs. 
The matrix C Q( g will then acquire a block-diagonal form, which facilitates the calculation of its inverse. Explicitly, 
the non-vanishing brackets are 

C\lx,2ky ■= [Xllx,X2ky]GP = *7feJ^(^ _ iJ) = C 2 kx,lly, 

C 3 x,4y : = [X3x, XiylcP = -AS(X ~y) = -C i3 ,Zy = -C 4 y,3x, (111.14) 
C5xfiy '■= [X5x, X6y\GP = A5(x — y) = —Cq$^ = —Cq^^x- 

The inverse matrices to these blocks are then 

G Ux,2ky = c 2kyMx = -^^(f _ fa 

CZ " AV ' = -^\hy^ +a ^ Ci3 '" V ' = ^\hj\ +a ^ (III.15) 

C 5x,6y = 1 + C 6x,5y = 1 + ^ 

4tt\x — y\ 47r|x — y\ 

Here ciij are completely arbitrary numbers. It may be surprising that the inverse matrices are not determined in 
a unique way and that they even do not inherit the symmetries of the matrices to which they are inverse. This is 
because these are infinite-dimensional matrices to which standard theorems of linear algebra do not apply. That the 
results (III. 15) are correct for arbitrary values of the constants can be verified by direct computations, e.g. 

J d 3 yC 3S ^C 4g . 3S =-J cPy A6(z- y) + a 3 4 J d?yA5{z- y) = -A, ( ^_ - ) = 8(2 - z). (111.16) 
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That the term proportional to 034 vanishes follows from the following reasoning: imagine that the area of integration 
is a bounded open subset of £1 C M 3 such that z € 0. Then the application of Gauss theorem (which is assumed to 
hold for distributions) allows to rewrite the integral of A5(z — y) as a flux of the gradient V^j5{z — y) through the 
boundary dQ. But the Dirac delta and its derivatives vanish everywhere except the points in which their argument is 
zero. Hence, the integration of VgS(z — y) with respect to y over the region that does not contain z (such as <9J1) must 
necessarily give zero. Now for any z € M 3 one can hnd Q, 3 z and decompose the integral over R 3 into the one over f2 
and ]R 3 \ f2. The integral over the complement of £1 vanishes, since there the argument of A<5 is always nonzero. This 
completes the argument. 

So, the arbitrary constants appeared in the inverse matrix C Q/3 . Will then the GDB also include them? The answer 
is no! Recall that the the purpose for the particular construction of GDB was to enable consistent replacement of 
the brackets by the commutators or anti-commutators of operators, dependently on whether the variables are even 
or odd. In order for this to be possible, the GDB was required to posses appropriate symmetries. This symmetries 
will not be satisfied if the inverse matrix (III. 15) is used in the construction, unless 

^33 = a 44 = »55 = a 66 = 0, 043 = — 1234, = — a$Q. (III. 17) 

So the freedom in the construction of the GDB is now reduced to the two parameters 034 and 056- The partially 
solved formula for the GDB that is convenient for performing farther calculations is then 

[F,G] GD = [F,G] GP + i-f kk > J d 3 x ([F,xikx]Gp[X2k'x,G} GP + [F,X2k'x]Gp[xikx,G} G p) 

+ ll d3xd3 y (44^ - fl34 ) {w^¥) [Gn2map - [F ^ mGP W^¥)) (IIL18) 

+ // d 3 xd 3 y( - ,i - +a 56 ) (^—{G,xe(y)} G p-lF,X6(y)}GP- ' 



Air\x-y\ J \SAq(x) 5A (x) 

However, another requirement for the GDB was that it is consistent with all the second class constraints. It appears 
that there is a problem with (III. 18) with this respect. To see this, take F — Aq(z) and G = J d 3 x\(x)ip(x), where A 
is arbitrary function that does not depend on the canonical fields. It follows from (III. 18) that 

[A Q (z),G] G D=ie Jd 3 x ( 4 ^ |a !_ |j + a 56 ) A(f)^(f). (111.19) 

But the constraint xe tells us that AA = ie(mp — ipn). which, under the physically reasonable assumption that A 
is bounded everywhere, can be integrated to yield 



t^^l^ + a^ ,111.20) 



where the constant doo does not contribute anything to the brackets and fact could be set to zero without any loss of 
generality in the subsequent analysis. If (III. 20) is used, the calculation of the bracket now gives 

\Mz), G] GD =iej i -|^A(f)V>(*)- (HI.21) 

Clearly, (III. 19) agrees with (III. 21) if and only if = 0! This is a strange result, since the bracket (III. 18) of any 
dynamical variable with any of the constraints Xi, i = 1,2, 3, 4, 5, 6 can be verified to vanish for ARBITRARY values 
of 034 and (J56. So it seams that it is consistent with the constraints! The point is that in obtaining (III. 20) I did not 
use the constraint xe but its integrated version 

Xe(x)d 3 x _ ^ f (nip- i(yjf) (x) 3 



My) ■■= / 2 V -1 = ie / I 1- -1 d x + ^ o( ^ ~ a °°- (IIL22) 

J 47r|x-y| J iTr\x-y\ 

The bracket (III. 18) of a dynamical variable with this integrated constraint does not vanish in general, unless as@ = 0. 
But how it can be that the bracket with xe vanishes but the bracket with the integral of x& does not? Does this 
fact contradict the linearity of the bracket? This apparent paradox can be traced back to the following calculation 
(certainly incorrect) 

'JiShJ^'-* - / ^/^tfqf- (,,L23) 
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The first equality stems from the fact that for any y £ M 3 we have J d 3 xA5(x — y) = (use Gauss theorem). That 
the result is wrong can be seen by substituting 



AS(x-y) 
An\z- y\ 



_d_ 

dy i 



1 



d 



Att\z — y\ dy'' 



S{y-x) - — - 



d_ 

dy l 



1 



4tt|z - y\ 



S(y 



d_d_ 

Qyl Qyl 



1 



4ir\z — y\ 



6(y-x) (HI.24) 



into the RHS of (III. 23). Then the integral of the first component of (III A) vanishes, since this component is the 
divergence of a vector held that vanishes everywhere except the point y = x. The integral of the second component 
of (HI A) gives —1, since this component is just equal to —S(z — y)S(y — x). on account of (III. 4). Hence, we finally 
obtain from (III. 23) that = —1, which seems incorrect. The mistake was made in the second equality of (III. 23), 
where I tacitly assumed that the order of integration with respect to x and y can be swapped. It appears that it can 
not. This is the reason why the linearity of the brackets that involve integrations in their structure cannot be naively 
exploited. 

In order to be able to use both the differential and integrated form of the constraints consistently with the bracket, 
I will set = = 0. Finally then the GDB does not contain any arbitrariness and is given by 

[F, G] GD = [F, G] GP + i 7 g fc , J d 3 x {[F, xiks} G p[X2k>s, G] GP + [F, X2k's]Gp[xikS, G]gp) + 
d 3 xd 3 y 



II 



d 



A-k\x — y\ \dy k 5-K k (y) 
or, even more explicitly, 

[F,G} GD ^[F,G} GP +z 7 ° kk , 
d SF d 8G 



SF 
[G,7 2 (x)] GP - [F, 72 (£)] G p 



SG 



J d 3 x([F, 



Qyk Sn k (y) SA (x) 



Xlkx]Gp[X2k'x, G] GP + [F, X2k'S\Gp[Xlkx, G] G p) 



F ^[G,x6(y)]GP - [F,x&(.y)]GP 



6Aq(x) 
(111.25) 



// 



d SF d SG 



SF 



dy k Sir k (y) dx 3 ' SAj(x) dxi SAj(x) dy k Sir k (y) SA (y) \Sir°(x 



SG 



+ie 



d SF 



^ + 



SF 



dy k Sir k {y) SA (y) 



[G, (nip - ipir)(x)] G p - ie[F, (irip - ipir)(x)} GP 



d 



SF 



d xd 3 y 
4tt\x — y\ 
SG 



Sir°(x) J SA (y) 



(111.26) 



SG 



^ + 



SG 



dy k SiT k {y) SA (y) 



From now on, the canonical variables A ai ir a will be referred to as electromagnetic variables, whereas ip, ip, tt and W 
as spinor variables. 

1) If both F and G depend on spinor variables only, then the only contribution to the bracket comes from the first 
line. The bracket is the same as the one calculated in the free Dirac field theory considered in [14]. 

2) Consider electromagnetic variables. It is clear that from (III. 26) that [A a (x), Ai,(y)] G D = [ir a (x), TT b (y)] G D = 0. 
The bracket of tt° with any other variable is certainly zero, since ir° is a constraint. The bracket of Ai with any 
variable G is 



{Mx),G] GD 



SG 



I i3 M 



d_ f SG \ 



Stt z (x) ' J ™ dx l \Ait\x — y\J dy k \Sir k (y)J ' 
from which the bracket for the canonically conjugate pair of electromagnetic variables follows 

1 



[^(x),7T J 



\GD 



Si5(x-y) + — — 



dx % dxi \Att\x — y\ 

The bracket of A with Gem that depends on the electromagnetic variables only is 



[A (y),G E M]GD = J d 3 x 



d 
dxi 



1 



d f SG 



EM 



4ir\x — z\dxi \Stt°(x) 



d 
dxl 



4ir\x — z\ J Sir°(x) 



SG 



EM 



(111.27) 



(111.28) 



(111.29) 



The integrand is a divergence of a vector field and hence the expression vanishes if SGem /Sir°(x) has bounded support, 
or just tends to zero sufficiently fast with x going to infinity. Note however that if, say, Gem = J \(x)ir° (x)d 3 x, where 
A does not tend to zero at infinity, then the bracket of G E m with A (y) will not vanish! So, the smeared constraints 
are consistent with the bracket (III. 26) if and only if the smearing functions decrease sufficiently rapidly at infinity. 
3) The bracket of A with Gsp that depends on the spinor variables only is 



[A {y),G S p] G D = ie 



d 3 3 



Att\x — y\ 



[G S p, {Trip - ipn)(x)] GP . 



(111.30) 



10 



4) The bracket of Fsp that depends on the spinor variables only with 7T- 7 (y) is 

[F S p,ir j (y)]GD = ie J d 3 x^- ( 47r |i_^ ) [F sp ,tt(x)^(x) ~^(x)W(x)] GP , (111.31) 

which is the same as —[Fsp,djA {y)]cD- This observation led Weinberg [3] to define a combined variable 

tt^ := n j + djA Q , (111.32) 

which has trivial bracket with spinor variables 

[F SP ,iri(if)}GD = 0. (111.33) 
Note that from X3 = and X6 = it follows that 

djiri = 0. (111.34) 
The bracket of 7r^_ with arbitrary dynamical variable F is 

SF f 3 ( d 5F \ ( d 1 



V^W™ = 6AM-J dX WsA^j) [dxJlW^l 1 ' (IIL35) 



from which it follows that 



lMx),<m GD - SfSix-y) + ^ . (111.36) 

B. The equations of motion 

Since the GDB is now used that is consistent with the second class constraints and all the constraints are now second 
class, we can freely use them to simplify the form of observables that are of interest. For example, the distinction 
between canonical, first class, total and extended Hamiltonians is now spurious: the simplest form of the Hamiltonian 
can be given as 

{Uem(x) + H D (x) + Ui{x)) d 3 x 7 

(111.37) 



H = J (H EM (x1+H D (i)+%(f))ft, 



The equations of motion for the canonical variables are 

My) = [My),H]GD = ni(y), (IH.38) 

¥ (y) = [k 1 (y) ,H] GD = d Fji (y) - e (^y ty) (y ) => & = V x B - J 

i>(y) = -^J^jipiy) - ieA (y)tp(y) - im-f°tp(y) {i^ a V a - m) ip = 

The equations certainly need to be supplemented by the constraints. The constraints \i an d X2 can be used just to 
eliminate the variables ir and W from the formalism once and for all. Similarly, Xb tells that ir° = is not a physical 
degree of freedom. The constraint X3 then reduces to the Gauss law Vn = etp^°ip and xe tells that Aq is not an 
independent variable but a functional of matter fields 

Mx)=e[$fpM d %. (111.39) 
J Air\x-y\ 

For simplicity, I imposed on Ao the condition of vanishing at infinity. The general bounded solution to the constraint 
is obviously (III. 20). The reader is encouraged to verify the effects produced by nonzero constant of (III. 20) in 
the final results. 

Note farther that (III. 32) allows for the elimination of tt^ in favor of 7r^ and ip, ip. Then 7r^_ can be eliminated in 
favor of Aj on account of the second equation of (III. 38). Finally, the only electromagnetic physically important fields 
are Aj's, which are restricted by the Coulomb gauge condition Xi = djA-* = 0, so we end up with the two degrees of 
freedom of the electromagnetic field, as desired. 
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C. Transition to the interaction picture 

In order to construct the quantum theory in the case of the free Dirac field [14], we had to find a general solution 
to the evolution equations for fields. The arbitrary operator coefficients in the general solution appeared to obey 
an exceptionally simple commutation rules with themselves and the Hamiltonian and other physically important 
observables such as momentum or electric charge operator (although these commutation relations were not explicitly 
verified in [14]). This simplicity was crucial and allowed the construction of Fock space carrying the representation of 
all the commutation relations following from the GDB. 

Even the first look on (III. 38) allows to see that no such simple solution to the system of equations for QED in 
the Coulomb gauge is possible. One could try to differentiate the second equation and then use the third in order 
to eliminate 7T* from the system. Also, the first equation can be used to eliminate Aq. This can be done, but the 
resulting relation between Aj and ip, ip is fairly too complicated to be exactly solvable. 

The way out of this difficulties is to pass to the interaction picture and calculate physical quantities perturbatively. 
To accomplish this, let us decompose the Hamiltonian (III. 37) into the free part Ho and the interaction part Hi 
according to 



H (t) = / CH E M{t,x) + H D (t,x))d 3 x, Hi(t) = / ni(t,x)d 3 x, H = H (t) + H T (t) 



(111.40) 



Up to this point, I did not even bother to write the time argument of fields. It was understood that ip{x), Aj(x), 
etc. depend on time and their evolution is determined by the total Hamiltonian. Now the Hamiltonian itself is time- 
independent, but its constituent parts, such as Ho and Hj do depend on time, as indicated in (III. 40). I shall now 
chose an instant of time, for simplicity let it be t = 0. From now on Ho and Hi should be understood as Ho(0) and 
Hj(0) (the same concerns Hp and Hem)- Let F(t) be a dynamical variable, whose time evolution is determined in a 
usual way, by H. By the interaction picture of F we will understand a variable F^(t) whose evolution is determined 
by H and whose value at t — is equal to the value of F at this instant: 



F ( i)(0) 
F ( i)(0) 



[F {I) (0),H ] GD , 

' F (im,Ho}_=[[F {I) (0),Ho] GD ,H 



GD 



' GD 



(111.41) 



^7(0) 



F, 



(n-l) 



(I) 



(0),H 



GD 



[[F {I) (0),H _ 



gd> H o\gd'" Ho 



GD 



From this definition it follows that 



F {I) (t)=[F {I) (t),H ] 



GD 



(111.42) 



for any t (use the Tylor expansion). Another immediate consequences of the definition are that the time derivative 
of an interaction picture variable also evolves in the interaction picture (the same is obviously true for spatial partial 
derivatives). Finally, for any variables F and G the following implication holds 



F (im = [ F (im,H ] GD = G(0) =► F {I) (t) = G (I) {t) 



(111.43) 



for any t. This latter corollary will be used repeatably below (the relations between the derivatives at t = will be 
obtained and the relations between the interaction picture fields at any time will then be assumed). 
Using (III. 37) and (III. 26), the GDB of any dynamical variable F with H can be calculated 



[F,H ]gd = [F,Ho\ G p + iliv j ([F,Xus}Gp[X2i>x,H D } G p + [F,X2i>x}Gp[xiis,H D } GP )d 3 x 
afxafy ( d SF 



+ ic 



SF 



4ir\x~y\ \dy k 5w k (y) 5A (y) 
This simplification follows from 

SHo 



(111.44) 



5H _ k 



dk lTT = °> 

0TT K 



SA 



[H D , (tt0 - ipn)(x)]GP- 



a i n ^0 r, r\ SH 



(111.45) 
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From (III. 44) it follows that if F$p depends on the spinor variables only then its bracket is the same as the one that 
would be obtained in the theory of the free Dirac field discussed in [14]. Hence, the interaction picture field 
satisfies the free Dirac equation 

{i^da - m) V>(/) = 0. (111.46) 

For other fields, the calculations give 

[A i (0,x),H ] GD = TT i L (0,x) => ii (J)i (i,£)=7ri (J) (i,£), 

[TT i 1 _{0,x),H ] GD = d j F ji (0,x) => ■k i L(I) {t,x)=d j F {I):ji {t,x), (HI 47) 



[A (0,x),H]cD = ~e j 47r | f _ ^ — d y A (i)o{t,x) = -e j 



A-k\x — y\ 



The implication (III. 43) was used. Note that the evolution of Aq in the interaction picture is nontrivial, contrary to 
what seems to have been suggested in [3]. However, the point is that A is no longer necessary, since the first and the 
second equations now give a simple equation for Arjy 

A (I) i=d j F (I)ji =AA (I)i => □i4 (J)i = 0. (HI-48) 

This equation, together with the Coulomb gauge constraint 

djA {I)j = (111.49) 

and the equation (III. 46) for are the only important equations for the evolution of interaction picture fields. Now 
these equations are sufficiently simple that the general solution can be readily found. Indeed, (III. 46) have already 
been solved in [14]. The general solution to (III. 48) and (III. 49) is 

A (I)i {x) = j dT k (e- ik - x e iX (k)b x (k) + e ihx e* x (k)b{(k)) , k'e^k) = 0. (111.50) 

The assumption is made that k° — \k\ and dTk = d 3 fc/27r 3 2fc°, in complete parallel with the conventions of [14], 
which were adapted to the case of massless particles. Hence, e\(k) for A = 1, —1 form a basis of the two-dimensional 
subspace of M 3 orthogonal to k and b\(k) are arbitrary coefficients. The coefficients b\(k) are classically the complex 
conjugates of b\(k), but the conjugation is denoted by f, since it will pass to Hermitian conjugation in the quantum 
theory, so that the operators ^4m, are self-adjoint. The letter b was used to denote the coefficients in order to 
distinguish them from the creation and annihilation operators for the Dirac field that were introduced already in [14]. 
A convenient choice of basis e\ is given by 



^ / cos 6 cos (j) =F i sin <j> \ / sin 8 cos <j) 

e±i(k) = — j= I cos#sin0± icos(f> , fc = |fc| sin 9 sin 4> | . (III. 51) 

— sin j \ cos 8 



V2 




It is also convenient to define the matrices 

2 / cos 8 cos (j) — i sin <f> cos 8 cos <j> + i sin <f> \ 1 I ± ± \ 

e(k) := (ei(fc), e 2 (fc) ) = —= cos^sin^ + icoscj) cos 8 sin <j> — i cos , e := —= \ i — i , (III. 52) 

V / V2 y _ s - m0 _ sin6 ) J V2 \ n n / 

which are related by the standard rotation 

(cos </> — sin (j) \ / cos 8 sin 8 

sin./) cos0 1 0). (111.53) 

1 / \ -sin0 cos0 

When these conventions are adopted, the relation (III. 50) can be explicitly inverted with respect to the coefficients 6, 
&t 



b{k) = J d^xe^elR-^k) (k°A {I) (x) + iA {I) (xf) , 
b\k) = J d 3 xe- ik - x elR- 1 (k) (k°A {I) (x) - iA {I) (x)y 



(111.54) 



Think of b(k) as a column vector with components b\(k) for A = 1,-1. Similarly, A^ should be thought of as a 
column with components ^4(/)i, A (i)2i A (i)3- 
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D. Quantization 



[F,G\ T =i[F,G] GD , (111.55) 

where the upper sign applies whenever at least one of the variables is even and the lower corresponds to the case of 
both variables being odd. The evolution of the Heisenberg and the interaction picture operators can now be written 
simply as 

F(t) = e ttS F(0)e- lt ", F {I) (t) = e ttS « F(0) e - ltS ° . (111.56) 
It follows that knowing the commutation relation at t = it is easy to obtain the one that is relevant for t ^ using 
[F(t),G(t)} T =e^[F(0),G(0)} T e-^, [F (I) (t), G (I) (t)] T = e itfi °[F (/) (0), G (I) (0)} T e-^<> . (111.57) 
Hence, the relation (III. 36) leads to 

[A {I)i (t,x),n j ±{I) (t,y)}_ = [A i (t,x),n{(t,fi}_=i6i6(x-fi+i^ i -^ ( 1 — ) . (111.58) 

Since A^^t, x) = n]_,j-.(t,x) on account of (III. 47), the relation (III. 58), together with (III. 54), can be used to derive 
the commutation relations of the operators corresponding to the coefficients b and defined by (III. 50) and (III. 51) 

[bx(k),b x >(%)]- = (2ir) 3 2k°6 xx ,6(k-k') (111.59) 

(note that I do not use hats b above the annihilation and creation operators). The physical interpretation of b and b^ 
can be inspected in exactly the same way as the meaning of a and a 1 * was established in [14]. Namely, the energy- 
momentum and spin density tensors should be constructed and their commutation relations with b and determined. 
It follows that k ought to be interpreted as the momentum and A as the helicity of the particle created by b\(k). 
Note, however, that this interpretation is valid in the asymptotic regions in which the interaction is sufficiently week 
that the time evolution can be satisfactorily approximated by the free Hamiltonian H n . This situation is perfectly 
relevant for the description of scattering experiments. Similarly, the operators a a (p), a c a (p) and their conjugates, 
whose commutation relations, obtained in [14], are 



a„(p),al(p')\ + = [a c a (p),a c J,(p')\ + = (2n) 3 2E p d^S(p- p'), 
[a„(p),a a ,(p')] + = K(p),a c a ,(p')] + = \a a (p),a c J,(p')] = a*, (?') 



= 0, 

+ 



(111.60) 



should be thought of as describing Dirac particles and their anti-particles in the asymptotic in and out regions. Of 
curse, the commutators of a's with &'s all vanish. This follows from the fact that both A{ and 7r^ = Ai have vanishing 
Dirac brackets with spinorial variables (see (III. 27) and (III. 35)). The Hilbert space carrying a representation of the 
commutation relations between the interaction picture fields, can now be simply constructed as the Fock space of a, 
a' , a c , a c \ b, w . The asymptotic in and out states will inhabit this space. The non-interacting vacuum (i.e. the 
lowest energy state for the free Hamiltonian Ho), from which this Fock space is constructed, satisfies 

a a (p)\0) = a%(p)\0) = b x (k)\0) - 0. (111.61) 

To proceed with this description, I shall use the well known formula for the perturbaive expansion of the S matrix 
elements (the so called Dyson series): 

00 (—i) N r r 

5(1,2,-..^1',2',...)= J2 Scrr / / d i x 1 ...d 4 x N {0\...a 2l a 1 ,T{:h(x 1 ):---:h(x N ):} a{ a\ . . . |0) 

N+0 ' •> J 

(111.62) 

(compare with (6.1.1) of [3]). Here 1,2, ••• denote the incoming particles, 1', 2', ••• are outgoing particles (these 
numbers are assumed to include the information about the momenta, spin projections, helicities and types of particles), 
T{} is the time ordering, h(x) is the interaction Hamiltonian density in the interaction picture and :: denotes normal 
ordering operation (i.e. all the annihilation operators occur on the right of the creation operators). Using (III. 37) 
and (III. 39), one obtains 

h(x) = Ha{x) + hc(x), 

/~ .~ \ e 2 f 5(x° - y") (— - \ /— - \ (III. 63) 

h A (x) = eA {I)j (x) ^ ( /)7^(/)J (x), h c (x) = — J d 4 y (V> ( j)7V(i)J (x) (V'(/)7°V'(/)J (y)- 
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IV. EXAMPLES 
A. Compton scattering 

The Compton scattering is the process in which a photon with initial four-momentum k and polarization 

e;(fc) = a+e itl (k) + a_e;_i(/c), |a_| 2 + |a+| 2 = 1 (IV.l) 

and the electron with initial four-momentum p and the projection of spin onto the third spatial axis a interact to 
produce the outgoing photon with four-momentum k' and polarization 

e , i (k') = a' + e i , 1 (k') + a'_e i ,- 1 (k'), |«'_| 2 + 1«^| 2 = 1 (IV.2) 

and the outgoing electron with four-momentum p' and the spin projection a'. The circular polarization of the photon 
corresponds to a_ = or a + = 0, whereas the linear polarization corresponds to |a_| = \a + \. 

Since the calculations are going to be lengthy, it is extremely useful to simplify the notation. I shall define 

/t :=«+&*(£) + a_6 t _ 1 (£), f:=ar + b 1 (k') + ar_b-i(k'), a* := 4(p), a' := a a ,(p'), 



4 +) := / dT pll e-^" x u Vall {f)a all {f), x ( f ] := / aT p „e^" x v la „ tf')a$, &'), *l ~ J+) ^ " (_) 



x 



(IV.3) 



4 +) : = / rfr fe » e - ife "^ iA (fc")&A(fc") ; 4 _) : = *i +)t > ^ : = 4 +) + 

M 3 .- yy. 

The interaction density in the interaction picture (III. 63) can now be rewritten as 
h{x) = h A (x) + h c (x), 

e 2 f S(x° - y°) (IV.4) 
h A (x) = eM 3 u , xj x] x v , h c (x) = — J d A y ^ 4 2$ 2/i' 

and the Dyson series (III. 62) for the Compton scattering is 

S(f,a^f',a')= ^vf" /■■■/ ^...^(Ol/VTI:^!):---:/!^):}^^^). (IV.5) 
The iV = term gives 

(Ol/VaV^O) = [a', a f ]+ [/',/*]_ = (2tt) 6 2E k 2E p 5(k - k')S(p- p")<W (a'* + a+ + a'*_aJ) . (IV.6) 

Assume now that we wish to calculate the probability amplitude of the scattering event in which the final momenta 
k', p' are at least slightly different than the initial momenta k and p. This assumption is certainly allowable, since it 
is up to us probability of which we wish to calculate! Under this assumption, the delta functions vanish and one gets 
no contribution from the N — term. In fact, since the wave packets that describe the incoming beams of particles 
are usually not ideally localized in the momentum space, the forward scattering term (IV.6) may contribute slightly 
to the measured values of the nontrivial scattering, but I will not discuss that kind of technical complications here. 
For N = 1 the two terms corresponding to Ha and he need to be evaluated 

i J d 4 x(0\fa' : h A (x) : O t/t|0) - i J d 4 x{0\fa'T{: h c (x) :} a t/+|0>. (IV.7) 

Note that he needs to be time-ordered due to its nonlocal character (the time ordering of Ha is not necessary) . It is 
easy to verify that the first component vanishes. Hence, the Coulomb component needs to be calculated. However, 
this term also appears to vanish, which follows from the fact that 

(0|a' :x\x lV \,y v : a+JO) = 0. (IV.8) 

For N = 2 the contribution is 

H) 2 



2! 



d 4 xd 4 y(0\fa'T{: h(x)::h(y) :} a t/t|0>. (IV.9) 
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But 

: h(x) : : h(y) := (: h A {x) : + : h c (x) :) (: h A (y) : + : h c (y) :) ^ 
= : /iaO) : : : + : : : h c (y) ■ + : /ic*(z) : : : + : ftc(i) = : ^c(y) : • 

From (IV. 4) it is clear that the first term is proportional to e 2 and the remaining terms are of higher order in e. If 
we wish to find the correction to the S matrix that is of lowest nontrivial order in the small coupling constant e, we 
should neglect all the terms but the first. A straightforward calculation then gives 

(0|/V : h A (x)::h A (y) : o t/+|0) = e 2 M' w M 3 nnl (0\f'a' : x iX \x v : : yj y\y n , : oV f |0) 

= e 2 M} v M 3 nn ,Wx iyj fW : x\x v : : y\y n , : o+|0) (IV.ll) 
= e 2 Mj v Mi n , + \f>, yj ][ Xi ,f1] + [f',f s [x^y^}) (0\a> : x\x v : : yly n , : o+|0) 

Here and below the bracket [, ] means the commutator for even variables and anti-commutator for odd ones (I drop 
the subscript =f). Up to now, I have pulled out all the electromagnetic operators in terms of commutators. The last 
factor contains the free vacuum expectation value of the fermionic operators and it will now be computed 

(0\a' : x\x v : : y\y n , : a f |0) 

= [a'.ar/lfy^.a+J^lari/yt |0> - [a'.ytl^.o+KOlarJy^lO) + [a',at]<0| : ^a:^ : : y+y„/ : |0) (IV.12) 
= [a',x\][y n ,,cJ)[x\+\y^] [a', y\][x v , <J][x™, y<7>] + [a', a^+\ yM][x™ , y^]. 

From (III. 59) and (III. 60) it follows that some of the brackets are proportional do the Dirac deltas between initial 
and final momenta 

[/',/+] - {2ir) 3 2E k 5(k - k') (a'* + a++ a'*_a_) , [a', at] = {2wf 2E p 6(p - $)5 aa ,. (IV.13) 

Recall that we are trying to calculate the probability amplitude of the scattering event in which the final momenta are 
different than the initial momenta (the discussion below the formula (IV.6)). Hence, all the terms that are proportional 
to [/',/^] or [a',a^] simply vanish. Then inserting (IV.12) into (IV.ll) and multiplying by —1/2 that was present in 
front of the integral in (IV. 9) yields 

(0|/V : h A (x) : : h A (y) : aV f |0) = ~ M tl' M L> 

([f\x i ][y j J^][a\xj}[y n ^a^\y^]-[f\x i }[ yj (IV.14) 

+ V\vj][xiji][a\xl]ly n ,,a1][xp\yM] [f\y j ][x i j1][a\yl][x l ,,a1][x\- ) \yP]) . 

In order to obtain the contribution to the S matrix, this expression needs to be time-ordered and integrated over 
d 4 xd 4 y. Let as perform this task for the first two terms of (IV.14). They can be rewritten as 

-y (\i n {x,y)[x\ + \y( + H] X nl (y, x)^, y^]) , X ln (x,y) := MUM^f 1 , Xi ][ yj , f^\a' ,x\,][y n , ,a)}. (IV.15) 
Time ordering and integration then yields 

D 2 



~ JJd i xd 4 y{9(x°-y") (\ ln (x, y)[x ( ^ , yWt] - \ nl (y , x)^ , y^}) 

+ 6(y° x°) (\ ln (y, x) [y^ , *(+>t] - \ nl ( x , y) [ y (">t , x (-)]\ \ 
rr f , ( IV - 16 ) 

= -e 2 Jj d A xd*y {8(x° y°)A ; „(x, y)[x\ + \ y^} %° x°)\ nl (x, y)]y^ , x^]} 
= ie 2 Jj d i xd i y\i n (x, y)A ln (x, y), 

where 9(t) is the step function taking value for t < 0, 1/2 for t = and 1 for t > 0, which can be conveniently 
represented as an integral 

9(t) = lim -— / ds, (IV.17) 

v ; e^o+ 2m / m s + is K ' 
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and 

A ln (x,y) = i (e(x° y°)[x\ + \y^} %° x°)[y^ , *<">]) (IV.18) 

is called the propagator for the Dirac field. Treating the last two terms of (IV. 14) in the same way, a similar expression 
is obtained 

ie 2 J J d i xd i yXm(x,y)A ln (x,y), \ ln (x,y) := M^M^lf , yj ][ Xi , f][a' ,x\,}[y n , ,o)] (IV.19) 

(note that A differs from A just by the interchange of Xi and yj). Hence, the total contribution to the S matrix that 
is proportional to e 2 is 



ie 2 



J I d 4 xd 4 y (\ ln (x, y) + \ ln {x, y)) A ln {x, y), (IV.20) 



with A and A given in (IV. 16) and (IV.19). The two terms of (IV.20) correspond to the Feynman diagrams FIG.l 
and FIG. 2. It is however important to see that the usage of Feynman diagrams and postulating Feynman rules is 
not really necessary, but only useful. In principle one could simply work out the commutation relations in the Dyson 
series, just as I did. This procedure certainly does not require the presence of Poincare invariance. 
Using (IV. 3), the Dirac propagator can be calculated and recast in a convenient form 



d 4 a N(n')e~ iq ^ x 'y^ 

A(x, y) = lim / - ^ 4 [q 2 h 2 . , N(q) = u(q)u\q) = (m+ fl 7 , { := q al a . (IV.21) 
£-5-0+ J y^Tt) m — q — is 



Hrere u(q) was assumed to be constructed according to the conventions adopted in [13], which are consequently used 
in this article. 

In order to calculate Aj„ and Aj„ explicitly, it is necessary to use 

[a',x\,\ = e^W). \Vn'o)] = e^vW, [f M = e ife 'V*(fc'), [ Vj ,p] = c^e^k), 
lf , ,y j ] = e ik 'ye'*(k'), [x u /t] = e^e^k). 

This results follow straightforwardly from (IV. 3). Inserting these results into the formulas for A and A given in (IV. 16) 
and (IV.19), one gets 

X ln (x,y) = e i(k'*-ky +P ' X - P y) [u{p')r{k')] all [7V(*)«(p)] n<T , 

\ ln (x,y) = e i(k'y-k x+P < x - P y) [u(p>)t(k)l n [ 7 V'*(fc')"b)]„ ff , ( IV - 23 ) 
f ■■= fe*. 

inserting (IV. 23) and (IV.21) into (IV.20) one obtains the e 2 contribution to the S matrix in the form 



(IV.24) 
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To arrive at this result, one should first perform the integrals w.r.t. d A x and d 4 y, which will produce delta functions 
from the exponents of (IV. 23). These delta functions can be used in the subsequent integration over d*q. One ends 
up with one delta function that simply expresses the four-momentum conservation. 

Note that the expressions in the denominators of (IV. 24), namely m 2 — (p + k) 2 = —2p-k and m 2 — (p— k') 2 = 2p-k', 
never vanish, because the plane that is perpendicular to a null vector does not contain any time-like vectors. Therefore, 
the infinitesimal parameter e can be simply set to zero. Another useful observation is that N(p + k)j = m+ fi+ ft 
and N(p~ k')j° = m+ p"— ft'. Using all of these, the e 2 contribution to the S matrix element for Compton scattering 
can be finally expressed as (2tt) 4: S 4 (p + k — p' — k')iM, where 

iM = -ie 2 u{p>){f{k') ^^ m+ m ^^r rw) «&>). av.25) 

This formula agrees with the one that is obtained by standard methods from Feynman rules. It has exactly the same 
form as the first formula of Chapter 5.5 of [8]. 



B. e + ,e — > p + , p scattering 

Another exemplary calculation will concern the probability amplitude for the production of a pair of muon and 
anti-muon from scattering of electron and positron. The modification of the Lagrangian that would allow for inclusion 
of many kinds of fermions is straightforward: any term of the form tpLtp, where L is a matrix-differential operator, 

ought to be replaced by ip^Lip( r \ with (r) labeling different kinds of particles. The two terms contributing to the 
interaction Hamiltonian are now 

h A (x) = h e A (x) + h» A (x), 

If dV(s° - y°) + + = + + I d*y8(x°-y°) 

2 J 47r|x-y| ° J Air\x - y\ 

p e (x) = e^ e (x) 7 > e (x), p»(x) = e^W/ffi). 

(IV.26) 

Let a\ a c * denote the creation operators of electron and positron and b, b c the annihilation operators of muon and 
anti-muon. These operators will have to be supplemented by additional indexes for the momenta and spin projections 
of the corresponding particles, but I shall skip this labels in the beginning. 
The N = 1 term in the Dyson series is 

-i J d i x(0\bb c T{: h A (x) : + : h c (x) :}a^a^\0). (IV.27) 

The first component including h A is equal to 

(0|66 c : h e A (x) : a^a^O) + (0\bb c : h^(x) : a c V|0) = 0. (IV.28) 

The vanishing of these two terms follows from the fact that in the first one can commute b c to the right, without 
producing any non-vanishing anti-commutators, whereas in the second it is possible to commute a cjf to the left. The 
Coulomb part of (IV.27) is, however, nontrivial and reads 

-ijf d 4 xd 4 y S -^^-(0\bb c T{: p e (x)p»(y) :}a ct a+|0). (IV.29) 

This terms constitutes the contribution of order e 2 following from the N = 1 term of the Dyson series. I shall return 
to this contribution later. 

The N — 2 term of the Dyson series is 

~JJ d 4 a;d 4 y(0|&6 c T{(: h A {x) : + : h c (x) :) (: h A (y) : + : h c (y) :)}a c V|0>, (IV.30) 

but the only term that is proportional to e 2 is 

-\fj d 4 xd A y(0\bb c T{: h A (x) :: h A (y) :}a c t a t|0). (IV.31) 
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FIG. 3: The diagram corresponding to (IV. 32). 

The remaining terms in (IV. 30) are of higher orders in e. The calculations that are similar to those performed in the 
case of Compton scattering can now be performed, which allow one to recast these expression as 



i J J d 4 xd i yA :jk (x,y)X jk (x,y), 

\*{x,y) := e 2 MiMUxl\a^[xt,,at][b,y^[b c ,yZ,), 
where again the short notation was used 

x? := ^ (e)t (x), xf, := $\x), y$ := ^%), j£, := 0(y), 
at:=4(p), a* :=<#(*?'), b:=b x (k), b c := b c x ,(k'), 



(IV.32) 



(IV.33) 



where now p, <r and p', a' are the momenta and spin projections of the ingoing electron and positron, whereas k, 
X and k', X' are the momenta and spin projections of the outgoing muon and anti-muon. The propagator for the 
electromagnetic field 

A jk (x,y) := i (0(x° y )^^} %° z^]) , (IV.34) 

(the definition of x\ +>> and x\~^ was already given in (IV. 3)) is equal to 



A jk (x,y) = - hm J — 4 , N M = 5 jk (IV.35) 



The expression (IV.32) corresponds to the Feynman diagram FIG. 3. If the calculations are continued, (IV.32) 
acquires a form 



(2^S(k + k '- P -p')^ r {6 jk V ^ [v{p'Wu{p)\ a , a [u(k)j k v(k')] xy , (iy 3g) 

q 2 = (k + k') 2 = (p + p') 2 , 

which differs from what would be obtained by the application of Lorentz-covariant Feynman rules to the diagram 
FIG. 3. It is now necessary to recall about the Coulomb term describing the e 2 contribution from N = 1 term of the 
Dyson series. This term is given by (IV. 29). Some simple manipulations employing the anti-commutation relations 
and invariance of the measure w.r.t. exchange of x and y allow to rewrite (IV. 29) as 

~ 2ie2 JJ ^\x%\ 5{x ° y0)d{x0 y Q )e- ix - {v+p,) e iy - (k+k,) W{p')u{p)] a , a [u\k)v{k')] xx , . (IV.37) 



2 



Using the identities 



in\x-y\ J (2tt) 3 |^| 2 ' in\x-y\ J (2tt) 4 |^*| 



Jq-(x-y) 



(IV.38) 
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and the fact that the step function obeys 8(0) = 1/2, one can finally rewrite the Coulomb contribution as 

(2n)W(k + k'-p- p')-pi- [vHp'Hp)] a , a [uHk)v(k')] xy . (IV.39) 
\k + k'\ z 

Adding this result to (IV. 36) and omitting (2n) 4 5 4: (k + k! — p — p') we get ultimately 

iM = ^r, ab W)l a u(p)] a , a [u(k) 7 b v(k')] xx , , q 2 = (k + k'f = {p + p'f. (IV.40) 
This final result is clearly covariant and agrees with other references, e.g. [8]. 

V. CONCLUSIONS 

Given a Lagrangian formulation of a field theory, there exists an algorithmic procedure for finding all the constraints, 
the Hamiltonian and the commutation relations of all the fields with respect to the GDB. If there are first class 
constraints present, one can try a gauge invariant method of quantization such as BRST quantization (which was not 
discussed here) or, alternatively, one can eliminate gauge freedom by imposing gauge conditions. A requirement of 
consistency of gauge conditions with time evolution needs to be imposed and it may lead to additional constraints. 
If all the gauge freedom is eliminated, the remaining constraints are all second class and should be incorporated in 
the construction of GDB. Then, in order to pass to quantum theory, one should in principle seek for representation 
of the final commutation relations of important fields with respect to GDB in a Hilbert space. 

All these steps were performed for the case of electrodynamics with fermions. The causal structure of space-time 
was not employed and, indeed, all these steps can still be performed if the electromagnetic interaction is replaced 
or supplemented by gravity. However, even in the case of electrodynamics, the resulting commutation relations 
appeared to be rather complicated and no obvious way of representing them in a Hilbert space was visible. In the 
simplest case of the free Dirac field, considered in the first article [14], it was possible to find explicit solution to the 
field equations. If the fields are constructed in such a way that they satisfy field equations automatically, then one 
does not need to bother about representing of their commutation relations with the Hamiltonian. What is more, the 
remaining commutation conditions between the functions that parametrize the exact solutions were sufficiently simple 
that the representation could be found for them. On the contrary, in the presence of electrodynamics the equations 
are to complicated to be solved explicitly. It is a way out of this problem in the case of electrodynamics, since the 
Hamiltonian decouples into the free and interaction part, the latter being proportional to a very small fine structure 
constant. This allows for the transition to the interaction picture, in which the equations and commutation relations 
are sufficiently simple. Then the perturbative quantization can be applied. In the case of gravity, no such obvious 
decoupling, which would lead to the simple interaction-picture equations, seems to be possible. On the other hand, 
the non-perturbative field equations are even more complicated then those of electrodynamics and hence the attempts 
to represent them, without any simplifications, in a Hilbert space seems hopeless. Certainly, one can decouple the 
metric tensor to the Minkowski part and the deviation, which is commonly done, and quantize only the deviation. 
Then the interaction picture can be obtained, but the background independence is sacrificed and the resulting theory 
is non-renormalizable. Hence, these problems are serious, and the lack of Poincare invariance in the presence of 
gravity is certainly not an important obstacle. 
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